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Abstract 

We review theoretical ideas, problems and implications of different models for neu- 
trino masses and mixing angles. We give a general discussion of schemes with three 
or more light neutrinos. Several specific examples are analyzed in some detail, par- 
ticularly those that can be embedded into grand unified theories. 
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1 Introduction 



There is by now convincing evidence, from the experimental study of atmospheric and 
solar neutrinos 0], for the existence of at least two distinct frequencies of neutrino 
oscillations. This in turn implies non vanishing neutrino masses and a mixing matrix, in 
analogy with the quark sector and the CKM matrix. So apriori the study of masses and 
mixings in the lepton sector should be considered at least as important as that in the 
quark sector. But actually there are a number of features that make neutrinos especially 
interesting. In fact the smallness of neutrino masses is probably related to the fact that u's 
are completely neutral (i.e. they carry no charge which is exactly conserved) and are Ma- 
jorana particles with masses inversely proportional to the large scale where lepton number 
(L) conservation is violated. Majorana masses can arise from the see-saw mechanism ||, 
in which case there is some relation with the Dirac masses, or from higher- dimensional 
non-renormalisable operators which come from a different sector of the lagrangian density 
than any other fermion mass terms. The relation with L non conservation and the fact 
that the observed neutrino oscillation frequencies are well compatible with a large scale 
for L non-conservation, points to a tantalizing connection with Grand Unified Theories 
(Cut's). So neutrino masses and mixings can represent a probe into the physics at GUT 
energy scales and offer a different perspective on the problem of flavour and the origin of 
fermion masses. There are also direct connections with important issues in astrophysics 
and cosmology as for example baryogenesis through leptogenesis M and the possibly non- 
negligible contribution of neutrinos to hot dark matter in the Universe. 

At present there are many alternative models of neutrino masses. This variety is 
mostly due to the considerable experimental ambiguities that still exist. The most crucial 
questions to be clarified by experiment are whether the LSND signal || will be confirmed 
or will be excluded and which solar neutrino solution will eventually be established. If 
LSND is right we probably need at least four light neutrinos, if not we can do with only the 
three known ones. Which solar solution is correct fixes the corresponding mass squared 
difference and the associated mixing angle. Another crucial unknown is the absolute scale 
of neutrino masses. This is in turn related to as diverse physical questions as the possible 
cosmological relevance of neutrinos as hot dark matter or the rate of neutrinoless double 
beta decay (0^/3/?). If neutrinos are an important fraction of the cosmological density, 
say Q u ~ 0.1, then the average neutrino mass must be considerably heavier than the 
splittings that are indicated by the observed atmospheric and solar oscillation frequencies. 
For example, for three light neutrinos, only models with almost degenerate neutrinos, with 
common mass \m u \ « 1 eV, are compatible with a large hot dark matter component, but 
in this case the existing bounds on Ovfifi decay represent an important constraint. On the 
contrary hierarchical three neutrino models (with both signs of Ami) have the largest 
neutrino mass fixed by \m\ ~ J Am^ m « 0.05 eV. In view of all these important questions 
still pending it is no wonder that many different theoretical avenues are open and have 
been explored in the vast literature on the subject. 

Here we will briefly summarize the main categories of neutrino mass models, discuss 
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their respective advantages and difficulties and give a number of examples. We illustrate 
how forthcoming experiments can discriminate among the various alternatives. We will 
devote a special attention to the most constrained set of models, those with only three 
widely split neutrinos, with masses dominated by the see-saw mechanism and inversely 
proportional to a large mass close to the grand unification scale M GUT . In this case one 
can aim at a comprehensive discussion in a GUT framework of all fermion masses. This 
is to some extent possible in models based on SU(5)x U(1)f or on SO(10) (we always 
consider SUSY GUT's). 



2 Neutrino Masses and Lepton Number Violation 

Neutrino oscillations imply neutrino masses which in turn demand either the existence 
of right-handed (RH) neutrinos (Dirac masses) or lepton number L violation (Majorana 
masses) or both. Given that neutrino masses are certainly extremely small, it is really 
difficult from the theory point of view to avoid the conclusion that L conservation must 
be violated. In fact, in terms of lepton number violation the smallness of neutrino masses 
can be explained as inversely proportional to the very large scale where L is violated, of 
order Mqut or even Mpi. 

Once we accept L non- conservation we gain an elegant explanation for the smallness 
of neutrino masses. If L is not conserved, even in the absence of heavy RH neutrinos, 
Majorana masses can be generated for neutrinos by dimension five operators of the 
form 

o 5 = (MMElk + h . c . , (i) 

with H being the ordinary Higgs doublet, U the SU(2) lepton doublets, A a matrix in flavour 
space and A a large scale of mass, of order Mqvt or Mp\. Neutrino masses generated by 
O5 are of the order m v « v 2 /A for A^ ~ 0(1), where v ~ 0(100 GeV) is the vacuum 
expectation value of the ordinary Higgs. 

We consider that the existence of RH neutrinos v c is quite plausible because all GUT 
groups larger than SU(5) require them. In particular the fact that v c completes the 
representation 16 of SO(10): 16=5+10+1, so that all fermions of each family are contained 
in a single representation of the unifying group, is too impressive not to be significant. At 
least as a classification group SO (10) must be of some relevance. Thus in the following 
we assume that there are both v c and L non-conservation. With these assumptions the 
see-saw mechanism || is possible. Also to fix notations we recall that in its simplest form 
it arises as follows. Consider the SU(3) x SU(2) x U(l) invariant Lagrangian giving rise 
to Dirac and v c Majorana masses (for the time being we consider the v Majorana mass 
terms as comparatively negligible): 

£ = -v cT y v (Hl) + ^v cT Mv c + h.c (2) 
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The Dirac mass matrix = y v v/y2, originating from electroweak symmetry breaking, 
is, in general, non-hermitian and non-symmetric, while the Majorana mass matrix M is 
symmetric, M = M T . We expect the eigenvalues of M to be of order Mqut or more 
because v c Majorana masses are SU(3)x SU(2)x U(l) invariant, hence unprotected and 
naturally of the order of the cutoff of the low-energy theory. Since all v c are very heavy 
we can integrate them away. For this purpose we write down the equations of motion for 
v c in the static limit, i.e. neglecting their kinetic terms: 

dC 

- — = y v {Hl) - Mv c = . (3) 

From this, by solving for z/ c , we obtain: 

v c = M- l y u {Hl) . (4) 

We now replace in the lagrangian, eq. (|2]), this expression for v c and we get the operator 
5 of eq. (0) with 

y = -yjM-V , (5) 

and the resulting neutrino mass matrix reads: 

m v = m^M _1 m£> . (6) 

This is the well known see-saw mechanism result [[|: the light neutrino masses are 
quadratic in the Dirac masses and inversely proportional to the large Majorana mass. 
If some v c are massless or light they would not be integrated away but simply added to 
the light neutrinos. Notice that the above results hold true for any number n of heavy 
neutral fermions R coupled to the 3 known neutrinos. In this more general case M is an n 
by n symmetric matrix and the coupling between heavy and light fields is described by the 
rectangular n by 3 matrix mo- Note that for m v rs \j Aml tm ~ 0.05 eV and m v w m 2 D /M 
with mu ~ v ~ 200 GeV we find M ~ 10 15 GeV which indeed is an impressive indication 
for Mqut- 

If additional non-renormalisable contributions to O5, eq. |l|, are comparatively non- 
negligible, they should simply be added. After elimination of the heavy right-handed 
fields, at the level of the effective low-energy theory, the two types of terms are equivalent. 
In particular they have identical transformation properties under a chiral change of basis in 
flavour space. The difference is, however, that in the see-saw mechanism, the Dirac matrix 
ran is presumably related to ordinary fermion masses because they are both generated by 
the Higgs mechanism and both must obey GUT-induced constraints. Thus if we assume 
the see-saw mechanism more constraints are implied. 



3 Baryogenesis via Leptogenesis from Heavy v c Decay 

In the Universe we observe an apparent excess of baryons over antibaryons. It is appealing 
that one can explain the observed baryon asymmetry by dynamical evolution starting from 
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an initial state of the Universe with zero baryon number (baryogenesis) . For baryogenesis 
one needs the three famous Sakharov conditions: B violation, CP violation and no thermal 
equilibrium. In the history of the Universe these necessary requirements can have occurred 
at different epochs. Note however that the asymmetry generated by one epoch could be 
erased at following epochs if not protected by some dynamical reason. In principle these 
conditions could be verified in the SM at the electroweak phase transition. B is violated 
by instantons when kT is of the order of the weak scale (but B-L is conserved), CP is 
violated by the CKM phase and sufficiently marked out-of- equilibrium conditions could 
be realized during the electroweak phase transition. So the conditions for baryogenesis at 
the weak scale in the SM superficially appear to be present. However, a more quantitative 
analysis |7j, shows that baryogenesis is not possible in the SM because there is not enough 
CP violation and the phase transition is not sufficiently strong first order, unless m# < 
80 GeV, which is by now completely excluded by LEP. In SUSY extensions of the SM, in 
particular in the MSSM, there are additional sources of CP violations and the bound on 
ran is modified by a sufficient amount by the presence of scalars with large couplings to 
the Higgs sector, typically the s-top. What is required is that m h ~ 80 — 110 GeV, a s-top 
not heavier than the top quark and, preferentially, a small tan (3. This possibility has by 
now become very marginal with the results of the LEP2 running. 

If baryogenesis at the weak scale is excluded by the data it can occur at or just below the 
GUT scale, after inflation. But only that part with |B — L| > would survive and not be 
erased at the weak scale by instanton effects. Thus baryogenesis at kT ~ 10 10 — 10 15 GeV 
needs B-L violation at some stage like for m v if neutrinos are Majorana particles. The 
two effects could be related if baryogenesis arises from leptogenesis then converted into 
baryogenesis by instantons [[|. Recent results on neutrino masses are compatible with 
this elegant possibility |§. Thus the case of baryogenesis through leptogenesis has been 
boosted by the recent results on neutrinos ||. 

In leptogenesis the departure from equilibrium is determined by the deviation from 
the average number density induced by the decay of the heavy neutrinos. The Yukawa 
interactions of the heavy Majorana neutrinos v c lead to the decays v c — > IH (with I a 
lepton) and v c — > IH with CP violation. The violation of L conservation arises from the 
AL = 2 terms that produce the Majorana mass terms. The rates of the various interaction 
processes involved are temperature dependent with different powers of T, so that the 
equilibrium densities and the temperatures of decoupling from equilibrium during the 
Universe expansion are different for different particles and interactions. The rates Tal^) 
of AL = 2 processes depend also on the neutrino masses and mixings, so that the observed 
values of the baryon asymmetry are related to neutrino processes. Precisely, Tal^) ~ 
T 3 /A 2 where A is the large scale that appears in eq. dTj) and also in the expression 
of light neutrino masses m u ~ v 2 /A. The out-of-equilibrium condition Tal^) < T exp , 
where T exp ~ T 2 /M Pl is the expansion rate of the Universe, leads to T ^ A 2 /M P i which 
then implies the relation (when correct proportionality factors and sum over flavours are 
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included) : 

i n 12 HpV 

E<^[0.2eV(^ 1 ^) 1 /Y (7) 

i 

What exactly is the temperature T which is relevant for leptogenesis depends on the 
thermal history of the early Universe and goes beyond the realm of neutrino physics. But 
if T ^ A 2 /Mpi and A ~ Mqut then the upper limit is significant and compatible with 
present neutrino data. 

If the RH neutrinos are thermally produced, then the mass of the RH neutrino that 
drives L violation is limited by the reheat temperature after inflation, which in turn is 
typically required not to exceed 10 8 — 10 10 GeV. This limit can be evaded if the RH 
neutrinos are instead produced by large inflaton oscillations during the preheating stage 
PH- 



4 Four (or More) Neutrino Models 



The LSND signal [|5| has not been confirmed by KARMEN [|Tl]. It will be soon double- 
checked by MiniBoone fll2"f. Perhaps it will fade away. But if an oscillation with Am 2 « 



1 eV 2 is confirmed then, in presence of three distinct frequencies for LSND, atmospheric 
[0] and solar neutrino oscillations, the simplest possibility is to introduce at least four 
light neutrinos. Since LEP has limited to three the number of "active" neutrinos (that is 
with weak interactions, or equivalently with non- vanishing weak isospin, the only possible 
gauge charge of neutrinos) the additional light neutrino(s) u s must be "sterile" , i.e. with 
vanishing weak isospin. Note that v c that appears in the see-saw mechanism, if it exists, 
is indeed a sterile neutrino, but a heavy one. 

A possibility to accommodate atmospheric, solar and LSND evidences for neutrino 
oscillations without introducing one or more sterile neutrinos is to invoke CPT violation 



13]. The required independent frequencies are provided by different neutrino and anti- 



neutrino masses and the fit to the present data has a good quality fH]j . 

A typical pattern of masses that works for 4-v models consists of two pairs of neutrinos 
] [T5| ] with mass separation between the two pairs, of order 1 eV, corresponding to the LSND 
frequency. The upper doublet is almost degenerate at m 2 of order 1 eV 2 being only split by 
(the mass squared difference corresponding to) the atmospheric (solar) v frequency, while 
the lower doublet is split by the solar (atmospheric) v frequency. An alternative to this 2-2 
spectrum is given by a 3-1 pattern with 1 being a nearly pure sterile neutrino separated 
by the LSND frequency from the 3. The 3-1 spectrum leads to a comparable (poor) 
overall quality of fit as the 2-2 pattern. These mass configurations can be compatible 
with an important fraction of hot dark matter in the universe. A complication is that 
the data appear to be incompatible with pure 1-v oscillations for v e — v s oscillations for 
solar neutrinos [T(J and for — v s oscillations for atmospheric neutrinos |T7[]. There are 



however (after SNO, marginally) viable alternatives. One possibility is obtained by using 
the large freedom allowed by the presence of 6 mixing angles in the most general 4-z/ 
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mixing matrix. If at least 4 angles are significantly different from zero, one can go beyond 
pure 2-u oscillations and, for example, for solar neutrino oscillations v e can transform into 
a mixture of v a and z/ s , where v a is an active neutrino, itself a superposition of z/ M and v T 
(mainly v T ) |T3|. A different alternative is to have many interfering sterile neutrinos: this 



is the case in the interesting class of models with large extra dimensions, where a whole 
tower of Kaluza-Klein neutrinos is introduced. This picture of sterile neutrinos from extra 
dimensions appears exciting and we now discuss it in some detail [[0| . 



The context is theories with large extra dimensions. Gravity propagates in all D 



dimensions (bulk), while SM particles live on a 4d brane. As well known | 20| , this can 
make the fundamental scale of gravity Mb much smaller than the Planck mass Mp\. In fact 
for D = 5 + 4 we have a geometrical factor V$, the volume of the compact dimensions, 
that suppresses gravity, so that 

(M D ) S V s = (M Pl /M D ) 2 , (8) 

and, as a result, can be as small as ~ 1 TeV. For neutrino phenomenology we need a 
really large extra dimension, with a radius R at least of the order of the scale set by the 
observed solar oscillation frequencies, 1/R ^ 0.01 eV or R ^ 0.02 mm. If we insist on 
having around 1 TeV we can assume, for instance, one compact dimension with radius 
R and 5 — 1 dimensions with a common radius R', such that the volume Vs = (27r) <5 RR' 5-1 
fits eq. (H). In string theories of gravity there are always scalar fields associated with 
gravity together with their SUSY fermionic partners (dilatini, modulini) pT|. These are 
particles that propagate in the bulk, have no gauge interactions and can well play the role 
of sterile neutrinos. The models based on this framework |22], ^] have some good features 
that make them very appealing at first sight. They provide a "physical" picture for v s . 
In the simplest case the theory includes a 5d fermion \l/(x, y), which decomposes into two 
4d Weyl spinors u s (x,y) and u' s (x,y) and contains a KK tower of recurrences of sterile 
neutrinos: 

!/.(*, y) = -/=E • (9) 

V ZnR n 

The tower mixes with the ordinary light active neutrinos in the lepton doublet I: 

S mix = [ d 4 x^=u s (x J 0)H(x)l(x) . (10) 



The interaction is restricted to the 4d brane at y — where SM fields live. Since the 5d 
spinor u s (x,y) has mass dimension 2, we need the mass parameter M 5 = M 5 d (2>kR') s - 1 to 
keep the Yukawa coupling constant dimensionless. From eqs. (§,[||IUp, after electroweak 
symmetry breaking, we find: 



Jd 4 xE^^\x)u a (x) , (11) 



where (H) = v/ a/2 and v a is the active neutrino embedded in /. Note that the geometrical 
factor M D /Mpi, which automatically suppresses the Yukawa coupling h, arises naturally 
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from the fact that the sterile neutrino tower lives in the bulk. An additional mass param- 
eter fi, related to a possible bulk mass term for the 5d fermion ty(x,y) is also allowed (in 
more realistic realizations, more 5d fields and L- violating interactions can be present). 

The pattern of oscillations results from the superposition of infinite components with 
increasing frequencies ~ n 2 and decreasing amplitudes ~ 1/n 2 . The leading oscillation fre- 
quency, v Am 2 , and the dominant mixing angle are determined by n and m = hvMp/Mpi, 
whereas the number of KK excitations that effectively take part in the oscillation is con- 
trolled by 1/R. Indeed, if 1/R 3> V Am 2 , the KK modes, whose masses are approximately 
given by n/R, decouple, with the possible exception of the lightest mode. If on the con- 
trary 1/R^, V Am 2 , then several KK levels participate to the oscillation and the resulting 
energy dependence of the survival/conversion probability can appreciably differ from that 
of the two level case. Indeed the contribution of a few KK states makes the solar oscilla- 
tion spectrum more compatible with the data. Note in passing that u s mixings must be 
small due to existing limits from weak processes, supernovae and nucleosynthesis p3|, so 



that the preferred solution for this KK v model is MSW SA. Instead the KK states should 
decouple in the case of atmospheric neutrino oscillations. These constraints fix the range 
of admissible values of R as specified above. 

In spite of its good properties there are problems with this picture, in our opinion. 
The first property that we do not like of models with large extra dimensions is that the 
connection with GUT's is lost. In particular the elegant explanation of the smallness of 
neutrino masses in terms of the large scale where the L conservation is violated in general 
evaporates. Since Md ~ 1 TeV is small, what forbids on the brane an operator of the 
form (Hl)J Xij(Hl) j / Mr> which would lead to by far too large v masses? One must impose 
by hand L conservation on the brane and that it is only broken by some Majorana masses 
of sterile z/'s in the bulk, which we find somewhat ad hoc. Another problem is that we 
would expect gravity to know nothing about flavour, but here we would need RH partners 



for z/ e , Vn and v T . Also a single large extra dimension has problems, because it implies [24 



a linear evolution of the gauge couplings with energy from 0.01 eV to ~ 1 TeV. But 
for more large extra dimensions the KK recurrences do not decouple fast enough. Perhaps 
a compromise at d=2 is possible. In conclusion the models with large extra dimension 
are interesting because they are speculative and fascinating but the more conventional 
framework still appears more plausible at closer inspection. 



5 Three-Neutrino Models 

We now assume that the LSND signal will not be confirmed, so that there are only two 
distinct neutrino oscillation frequencies, the atmospheric and the solar frequencies. These 
two can be reproduced with the known three light neutrino species (for other reviews of 
three neutrino models see 

Neutrino oscillations are due to a misalignment between the flavour basis, v' = (u e , z/ M , u T ), 
where v e is the partner of the mass and flavour eigenstate e~ in a left-handed (LH) weak 



H). 
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isospin SU(2) doublet (similarly for and u T )) and the mass eigenstates v = (^1,^2,^3) 

v' = Uu , (12) 

where U is the unitary 3 by 3 mixing matrix. Given the definition of U and the transfor- 
mation properties of the effective light neutrino mass matrix m u : 

v lT m v v' = v T U T m v Uv (13) 
U T m u U = Diag (777,1,777,2,777,3) = m diag , 

we obtain the general form of m u (i.e. of the light v mass matrix in the basis where the 
charged lepton mass is a diagonal matrix): 

m v = Um diag U T . (14) 

The matrix U can be parameterized in terms of three mixing angles #12, #23 and #13 
(0 < 9ij < 7r/2) and one phase <p (0 < ip < 2ti) P5| , exactly as for the quark mixing 
matrix Vckm- The following definition of mixing angles can be adopted: 



U 




up 



f C13 s 13 e 

1 I -, l2 r 12 I (15) 

V-sise"^ da 




where Sjj = sin%, q,,- = cos In addition we have the relative phases among the 
Majorana masses mi, m2 and m^. If we choose 777,3 real and positive, these phases are 
carried by m 12 = |m 12 |e i<?il ' 2 . Thus, in general, 9 parameters are added to the SM when 
non vanishing neutrino masses are included: 3 eigenvalues, 3 mixing angles and 3 CP 
violating phases. 

In our notation the two frequencies, Amj/AE (I = sun,atm), are parametrized in 
terms of the v mass eigenvalues by 

Am Ln = l Am L| ; Am ltm = \ Am ls\ ■ (16) 

where Am^ = | r7T, 2 1 2 — |^i| 2 and Am^ = m| — |m 2 | 2 . The numbering 1,2,3 corresponds 
to our definition of the frequencies and in principle may not coincide with the ordering 
from the lightest to the heaviest state. From experiment, see table [T], we know that 

C23 ~ S23 ~ 1/ y/2, corresponding to nearly maximal atmospheric neutrino mixing, and 
that S13 is small, according to CHOOZ, s 13 < 0.2 [p9fl . The solar angle 9 12 is probably 
large (MSW LA, LOW, VO solutions) or even maximal for LOW and VO, but could, 
alternatively, be very small s\ 2 ~ O(10~ 3 ) ||16|| , if the now disfavoured MSW SA solution 
is also kept in our list. If we take maximal s 2 3 and keep only linear terms in u = s\^e lip 
from experiment we find the following structure of the Ufi (/ = e,/i,r, % = 1,2,3) mixing 
matrix, apart from sign convention redefinitions: 



U 



fi 



I Ci2 S X2 U 

-(si2 + d 2 n*)/v / 2 (d2 - s l2 u*)/y/2 l/y/2\ (17) 
V (S12 - c 12 u*)/V2 -(c 12 + s 12 u*)/V2 l/y/2, 
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Given the observed frequencies and our notation in eq. fllHD, there are three possible 
patterns of mass eigenvalues: 



Degenerate 
Inverted hierarchy 
Normal hierarchy 



1 771-1 1 ~ |m-2 1 ~ I ^3 1 3> — TTlj\ 

\mi\ ~ | TTi2 1 |rri3| 
1*713 | > \m 2 ,i\ 



In the following we will discuss the phenomenology for these different cases and the re- 
spective advantages and problems. 

5.1 Degenerate Neutrinos 

For degenerate neutrinos the average m 2 is much larger than the splittings. At first sight 
the degenerate case is the most appealing: the observation of nearly maximal atmospheric 



Table 1: Square mass differences and mixing angles |16|, [17], 18 . 





lower limit 

(3a) 


best value 


upper limit 

(3a) 


(A<Jla ( 10 ~ 5 eV2 ) 


2.3 


5 


37 


(AmLjLow (10- 8 eV 2 ) 


3.5 


8 


12 


A»L (10- 3 eV 2 ) 


1 


3 


6 


(tan 2 ) la 


0.24 


0.4 


0.89 


(tan 2 #12) low 


0.43 


0.6 


0.86 


tan 2 #23 


0.33 


0.8 


3.3 


tan 2 6>i3 








0.07 
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neutrino mixing and the experimental indication that also the solar mixing is large (at 
present the MSW SA solution of the solar neutrino oscillations appears disfavoured by 
the data ||16|| ) suggests that all v masses are nearly degenerate. Moreover, the common 
value of \m v \ could be compatible with a large fraction of hot dark matter in the universe 
if \m u \ ~ 1 — 2 eV. In this case, however, the existing limits [|3D| on the absence of (V/3/3 
(|m ee | < 0.2 eV or to be more conservative \m ee \ < 0.3 — 0.5 eV) imply []3l| double maximal 
mixing (bimixing) for solar and atmospheric neutrinos. In fact the quantity which is bound 
by experiments is the 11 entry of the v mass matrix, which in general, from eqs. ( |13"D and 



(|15"D, is given by : 



\m e 



sj 3 ) {m x c\ 2 + m 2 s\ 2 ) + m 3 e 



2i<t> „2 I 
s 13 1 



(19) 



which in this particular case (m^ cannot compensate for the smallness of s 2 3 ) approximately 
becomes: 



\m e 



\micl 2 + m 2 s\ 2 \ ^ 0.3 — 0.5 eV 



(20) 



To satisfy this constraint one needs m\ pa —m 2 (recall that a relative phase <p 2 — 0i is 
allowed between m x and m 2 ) and c 2 2 ~ s\ 2 to a good accuracy (in fact we need sin 2 29 12 > 



0.96 in order that |cos2# 
following texture 



12| 



cos 2 #12 — sin 2 # 12 | < 0.2). This is exemplified by the 



/ 



m,. 



m 



■l/y/2 l/y/2 \ 



-l/y/2 (l + 77)/2 (l + r/)/2 
V l/y/2 (l + 77)/2 (l + v)/2 J 



(21) 



where r\ <C 1, corresponding to an exact bimaximal mixing, S13 = and the eigenvalues 
are m,\ = m, m 2 = —m and m 3 = (l + rj)m. This texture has been proposed in the context 
of a spontaneously broken SO (3) flavor symmetry and it has been studied to analyze the 
stability of the degenerate spectrum against radiative corrections | 3"2|| . A more realistic 
mass matrix can be obtained by adding small perturbations to m v in eq. (2T): 

/ S 



rn 



-1/V2 
\ 0-e)/y/2 



-1/V2 
(l + 77)/2 
[l + v -e)/2 



(l-e)/V2 \ 
(l + ri -e)/2 
[l + v -2e)/2 J 



(22) 



where e parametrizes the leading flavor- dependent radiative corrections (mainly induced 
by the r Yukawa coupling) and 5 controls m ee . Consider first the case 5 <C 6. To first 
approximation # 12 remains maximal. We get Am 2 aun pa m 2 e 2 /r) and 







( Am 



13 



2 1 
sun 



1/2 



) 2 

L atm , 



\Am\ 

If we instead assume 5 ^> e, we find Am 2 



m ee <C m 



m 4 



1/2 



(23) 



2m 2 S, #23 ~ tt/4, sin 2 2#i2 pa 1 - S 2 /A. Also 



in this case the solar mixing angle remains close to tt/4. We get: 



'13 



2m 



(24) 
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too small for detection if the average neutrino mass m is around the eV scale. This 
example shows that there is no guarantee for m ee to be close to the range of experimental 
interest, even with degenerate neutrinos where the involved masses are much larger than 
the oscillation frequencies. However, an almost maximal solar mixing angle such as the one 
implied by the previous analysis, is difficult to reconcile with the MSW LA solution. Of 
course the strong constraint s\ 2 = c i2 can be relaxed if the common mass is below the hot 
dark matter maximum. It is true in any case that a signal of Qv(3(3 near the present limit 
(like a large relic density of hot dark matter) would be an indication for nearly degenerate 
z/s. 

In general, for naturalness reasons, the splittings cannot be too small with respect to 
the common mass, unless there is a protective symmetry |32|, [33|]. This is because the 
wide mass differences of fermion masses, in particular charged lepton masses, would tend 
to create neutrino mass splittings via renormalization group running effects even starting 
from degenerate masses at a large scale. For example, for m ~ 1 eV, the VO solution for 
solar neutrino oscillations would imply Am/m ~ 10~ 9 — 10~ n which is difficult to obtain. 
Even in the previous example, where, for 5 <C e, the corrections to Am 2 sun are quadratic in e 
rather than linear, we would need e < (10~ 3 /m(eV)) 2 in order to have Am 2 sun < 10~ 9 eV 2 . 
In this respect the MSW LA or LOW solutions would be favoured, but, if we insist that 
\m v \ ~ 1 — 2 eV, it is not clear that the mixing angle preferred by the data is sufficiently 
maximal. Summarizing, degenerate models with |m| ~ 1 — 2 eV as required if u's are a 
cosmologically important source of hot dark matter have some problems related to OuP/3 
limits and to naturalness. In comparison degenerate models with sub-eV common mass 
appear simpler to realize. 

It is clear that in the degenerate case the most likely origin of v masses is from 
some dimension 5 operators (Hl)J Xij(Hl)j / A not related to the see-saw mechanism m v = 
m^M -1 m£>. In fact we expect the v Dirac mass not to be degenerate like for all other 
fermions and a conspiracy to reinstate a nearly perfect degeneracy between mjj and M, 
which arise from completely different physics, looks very unplausible (see, however, [34|). 
Thus in degenerate models, in general, there is no direct relation with Dirac masses of 
quarks and leptons and the possibility of a simultaneous description of all fermion masses 



within a grand unified theory is more remote [35 



The degeneracy of neutrinos should be guaranteed by some slightly broken symmetry. 
Models based on discrete or continuous symmetries have been proposed. For example 
in the models of ref. |36[ the symmetry is SO (3). In the unbroken limit neutrinos are 
degenerate and charged leptons are massless. When the symmetry is broken the charged 
lepton masses are much larger than neutrino splittings because the former are first order 
while the latter are second order in the electroweak symmetry breaking. 

A model which is simple to describe but difficult to derive in a natural way is one 
37p where up quarks, down quarks and charged leptons have "democratic" mass matrices, 
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with all entries equal (in first approximation): 



nit 



(I 
1 

u 



+ Srrij 



(25) 



where rhj (/ = u,d, e) are three overall mass parameters and 5m f denote small perturba- 
tions. If we neglect Srrif, the eigenvalues of are given by (0, 0, 3 rhf). The mass matrix 
m/ is diagonalized by a unitary matrix Uf which is in part determined by the small term 
5m f. If 5m u ~ 5ma, the CKM matrix, given by Vqkm = U^Ud, is nearly diagonal, due 
to a compensation between the large mixings contained in U u and When the small 
terms 5mj are diagonal and of the form 8m f = Diag(— e/, e/, 5f) the matrices £// are 
approximately given by (note the analogy with the quark model eigenvalues 7r°, tj and t/): 



U] 



(1/V2 
1/V6 
VVV3 



-1/V2 
1/V6 
1/^3 




(26) 



At the same time, the lightest quarks and charged leptons acquire a non-vanishing mass. 



The leading part of the mass matrix in eq. (25) is invariant under a discrete S^l x S^r 
permutation symmetry. The same requirement leads to the general neutrino mass matrix: 
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1 
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m u = m 
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_l_ r 
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1 








V 








u 
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1 ) 





+ 5m v 



D2 



(27) 



where 5m v is a small symmetry breaking term and the two independent invariants are 
allowed by the Majorana nature of the light neutrinos. If r vanishes the neutrinos are 
almost degenerate. In the presence of 5m v the permutation symmetry is broken and the 
degeneracy is removed. If, for example, we choose 5m u = Diag(0, e, 77) with e < r\ <^ 1 
and r<e, the solar and the atmospheric oscillation frequencies are determined by e and 
7], respectively. The mixing angles are almost entirely due to the charged lepton sector. 
A diagonal 5m e will lead to a neutrino mixing matrix U ~U\ characterized by an almost 
maximal #12 , tan 2 #23 ~ 1/2 and 



e 



13 



\fmjm p 



(28) 



By going to the basis where the charged leptons are diagonal, we can see that m ee is close 
to m and independent from the parameters that characterize the oscillation phenomena. 

The parameter r receives radiative corrections that, at leading order, are logarith- 
mic and proportional to the square of the r lepton Yukawa coupling. It is important to 
guarantee that this correction does not spoil the relation r C e, whose violation would 
lead to a completely different mixing pattern. This raises a 'naturalness' problem for the 
LOW and VO solutions. We conclude by stressing that a non-vanishing Am atm , maximal 
612, large 623 and vanishing # 13 are not determined by the symmetric limit, but only by a 
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specific choice of the parameter r and of the perturbations that cannot be easily justified 
on theoretical grounds. It would be desirable to provide a more sound basis for the choice 
of the small terms in this scenario that is quite favourable to signals both in the 0u2/3 
decay and in sub-leading oscillations controlled by 6>i 3 . 

Anarchical models can be considered as particular cases of degenerate models with 



m ~ Am atm . In this class of models one assumes that all mass matrices are structureless 
in the leptonic sector. At present the data appear to indicate the MSW LA solution as 
the most likely one. For this solution the ratio of the solar and atmospheric frequencies is 
not so small: typically (Am 2 un ) LA / Am atm ~ 0.01 — 0.2 and two out of three mixing angles 
are large. One important observation is that the see-saw mechanism tends to enhance the 
ratio of eigenvalues: it is quadratic in mc so that a hierarchy factor / in becomes 
f 2 in m u and the presence of the Majorana matrix M results in a further widening of 
the distribution. Another squaring takes place in going from the masses to the oscillation 
frequencies which are quadratic. As a result a random generation of the vhd and M 
matrix elements leads to a distribution of (Am^ n ) M /Am^ (m that peaks around 0.1. At 
the same time the distribution of sin 2 9^ is peaked around 1 for all three mixing angles. 
Clearly the smallness of #13 is problematic. This can be turned into the prediction that in 
anarchical models #13 must be near the present bound (after all the value 0.2 for sin# 13 
is not that smaller than the maximal value 0.701). In conclusion there is a non negligible 
probability that if the MSW LA solution is realized and #13 is near the present bound than 
the neutrino masses and mixings, interpreted by the see-saw mechanism, just arise from 
structureless underlying Dirac and Majorana matrices. 

5.2 Inverted Hierarchy 

The inverted hierarchy configuration \ra%\ ~ \m,2\ 3> 1 777.3 1 consists of two levels mi and 1712 
with small splitting Am 2 2 = Am 2 un and a common mass given by \m\ 2 \ ~ |Am 2 tm | ~ 
3- 10~ 3 eV 2 (no large hot dark matter component in this case). One particularly interesting 



example of this sort ||40|| , which leads to double maximal mixing, is obtained with the phase 



choice m% = —m 2 so that, approximately: 

m diag = Diag(v / 2m, -y/2m, 0) . (29) 

The effective light neutrino mass matrix 

m v = Um diag U T , (30) 

which corresponds to the mixing matrix of double maximal mixing C12 = S12 = V v^2 and 
S13 = u = in eq. ([IT]). 

1/V2 l/y/2 \ 
U fi = I -1/2 1/2 l/y/2 , (31) 
1/2 -1/2 l/y/2j 



13 



is given by: 



m 








(32) 



The structure of m u can be reproduced by imposing a flavour symmetry L e — — L T 
starting either from (Hl)f\ij(Hl)j/A or from RH neutrinos via the see-saw mechanism. 
The 1 — 2 degeneracy remains stable under radiative corrections. The preferred solar 
solutions are VO or the LOW solution. The MSW LA could be also compatible if the 
mixing angle is large enough. 

The leading texture in fl32|) can be perturbed by adding small terms: 



m 



( 5 



\ 



1 



1 \ 

7] 7] 
V V J 



(33) 



where 5 and r\ are small (<C 1), real parameters defined up to coefficients of order 1 that can 
differ in the various matrix elements. The perturbations leave Am 2 4rrt and #23 unchanged, 
in first approximation. We obtain tan 2 #i 2 ~ 1 + 5 + r) and Am 2 nn /Am 2 tm « rj + S, where 



coefficients of order one have been neglected. Moreover #13 « 77. If rj 3> 5, we have 



'13 



Am2 sun 



< \ Ami, 



1 

Am 2 \ 2 



A m atm , 



(34) 



In the other case, ijCiwe obtain: 



a „ Am2 sun 



7 13 



Am 2 



a t in 



(35) 



There is a well-known difficulty of this scenario to fit the MSW LA solution |^0|, f|l]. 
Indeed, barring cancellation between the perturbations, in order to obtain a Am 2 un close 
to the best fit MSW LA value, 77 and 5 should be smaller than about 0.1 and this keeps 
the value of sin 2 2#i2 very close to 1, somewhat in disagreement with global fits of solar 
data |TJJ. Even by allowing for a Am 2 Mn in the upper range of the MSW LA solution, or 
some fine-tuning between rj and 5, we would need large values of the perturbations to fit 
the MSW LA solution. On the contrary, the LOW solution can be accommodated, but, 
in this case, # 13 and m ee estimated in ( |3~4"1 , |3~5D are too small to be detected by planned 
experiments. 



5.3 Normal Hierarchy 

We now discuss the class of models which we consider of particular interest being the 
most constrained framework which allows a comprehensive combined study of all fermion 
masses in GUT's. We assume three widely split z/'s and the existence of a RH neutrino 
for each generation, as required to complete a 16 dimensional representation of SO (10) for 
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each generation. We then assume dominance of the see-saw mechanism m u = m T D M~ x mD. 
We know that the third-generation eigenvalue of the Dirac mass matrices of up and down 
quarks and of charged leptons is systematically the largest one. It is natural to imagine 
that this property could also be true for the Dirac mass of z/'s: rn^ 9 ~ Diag(0, 0, mm)- 
After see-saw we expect m v to be even more hierarchical being quadratic in mo (bar- 
ring fine-tuned compensations between mu and M). The amount of hierarchy, m\jm\ = 
Am 2 ,^/ Am 2 un , depends on which solar neutrino solution is adopted: the hierarchy is max- 
imal for VO and LOW solutions, is moderate for MSW in general and could become quite 
mild for the upper Am 2 un domain of the MSW LA solution. A possible difficulty is that 
one is used to expect that large splittings correspond to small mixings because normally 
only close-by states are strongly mixed. In a 2 by 2 matrix context the requirement of 
large splitting and large mixings leads to a condition of vanishing determinant and large 
off-diagonal elements. For example the matrix 



x 2 



x 



(36) 



has eigenvalues and 1 + x 2 and for x of 0(1) the mixing is large. Thus in the limit 
of neglecting small mass terms of order 7711,2 the demands of large atmospheric neutrino 
mixing and dominance of 777,3 translate into the condition that the 2 by 2 sub determinant 
23 of the 3 by 3 mixing matrix approximately vanishes. The problem is to show that 
this vanishing can be arranged in a natural way without fine tuning. Once near maximal 
atmospheric neutrino mixing is reproduced the solar neutrino mixing can be arranged to 
be either small or large without difficulty by implementing suitable relations among the 
small mass terms. 

It is not difficult to imagine mechanisms that naturally lead to the approximate van- 
ishing of the 23 sub-determinant. For example |42|, ^3[ assumes that one v c is particularly 
light and coupled to /i and r. In a 2 by 2 simplified context if we have 

. ^-(?!) ■ — • 

then for a generic mo we find 

t , , 1 1 ( a 2 ab\ . . 

= mlM-'mn ^ - j^ ab fe2 J . (38) 

A different possibility that we find attractive is that, in the limit of neglecting terms of 
order m 12 and, in the basis where charged leptons are diagonal, the Dirac matrix 77773, 
defined by v c mr>v, takes the approximate form, called "lopsided" |0, [0| |46| : 





777 D oc . (39) 
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This matrix has the property that for a generic Majorana matrix M one finds: 




m T D M- x m D oc x 2 x\ . (40) 



The only condition on M _1 is that the 33 entry is non zero. However, when the approxi- 
mately vanishing matrix elements are replaced by small terms, one must also assume that 
no new 0(1) terms are generated in m v by a compensation between small terms in mo and 
large terms in M _1 . It is important for the following discussion to observe that given 
by eq. ( |39"D under a change of basis transforms asm^-* V^m^U where V and U rotate 



the right and left fields respectively. It is easy to check that in order to make mo diagonal 
we need large left mixings (i.e. large off diagonal terms in the matrix that rotates LH 
fields). Thus the question is how to reconcile large LH mixings in the leptonic sector with 
the observed near diagonal form of Vckm, the quark mixing matrix. Strictly speaking, 
since Vckm = U^Ud, the individual matrices U u and Ud need not be near diagonal, but 
Vckm does, while the analogue for leptons apparently cannot be near diagonal. However 
for quarks nothing forbids that, in the basis where m u is diagonal, the d quark matrix has 
large non diagonal terms that can be rotated away by a pure RH rotation. We suggest 
that this is so and that in some way RH mixings for quarks correspond to LH mixings for 
leptons. 

In the context of (SUSY) SU(5) there is a very attractive hint of how the present 
mechanism can be realized E7], |48j. In the 5 of SU(5) the d c singlet appears together with 
the lepton doublet (v, e). The (u, d) doublet and e c belong to the 10 and v c to the 1 and 
similarly for the other families. As a consequence, in the simplest model with mass terms 
arising from only Higgs pentaplets, the Dirac matrix of down quarks is the transpose of 
the charged lepton matrix: = [mi) T . Thus, indeed, a large mixing for RH down quarks 
corresponds to a large LH mixing for charged leptons. At leading order we may have the 
lopsided texture: 

/0 0\ 

m d = {nnf = 1 \v d . (41) 

\o o 1) 

In the same simplest approximation with 5 or 5 Higgs, the up quark mass matrix is 
symmetric, so that left and right mixing matrices are equal in this case. Then small 
mixings for up quarks and small LH mixings for down quarks are sufficient to guarantee 
small Vckm mixing angles even for large d quark RH mixings. It is well known that a 
model where the down and the charged lepton matrices are exactly the transpose of one 
another cannot be exactly true because of the e/d and fi/s mass ratios. It is also known 
that one remedy to this problem is to add some Higgs component in the 45 representation 
of SU(5) PSfl . But the symmetry under transposition can still be a good guideline if we are 
only interested in the order of magnitude of the matrix entries and not in their exact values. 
Similarly, the Dirac neutrino mass matrix is the same as the up quark mass matrix in 
the very crude model where the Higgs pentaplets come from a pure 10 representation of 
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SO(10): tud = m u . For the dominance of the third family eigenvalue as well as a near 
diagonal form could be an order of magnitude remnant of this broken symmetry. Thus, 
neglecting small terms, the neutrino Dirac matrix in the basis where charged leptons are 
diagonal could be directly obtained in the form of eq. fl3"9|). 

To get a realistic mass matrix, we allow for deviations from the symmetric limit in 
(f40|), where we take x — 1. For instance, we can consider those models where the neutrino 
mass matrix elements are dominated, via the see-saw mechanism, by the exchange of 
two right-handed neutrinos [43|. Since the exchange of a single RH neutrino gives a 



successful zeroth order texture, we are encouraged to continue along this line. Thus, we 
add a sub-dominant contribution of a second RH neutrino, assuming that the third one 
gives a negligible contribution to the neutrino mass matrix, because it has much smaller 
Yukawa couplings or is much heavier than the first two. The Lagrangian that describes this 
plausible subset of see-saw models, written in the mass eigenstate basis of RH neutrinos 
and charged leptons, is 



C 



M' 



-v 



c/2 



leading to 



/ \ ViVj . ViVj 

[m v )j~ oc — - H - 

V uh3 M M' 



(42) 



(43) 



where i, j = {e, [i, r}. In particular, if y e <C y^ y T and y' ~ y' T , we obtain: 



(6 



m,. 



m 



e 1 + rj 1 + r] 
e 1+7] 1+7] ) 



(N) , 



(44) 



where coefficients of order one multiplying the small quantities 5, e and r\ have been 
omitted. The mass matrix in fl44D does not describe the most general perturbation of the 
zeroth order texture (|40D. We have implicitly assumed a symmetry between and v T 
which is preserved by the perturbations, at least at the level of the order of magnitudes. 
The perturbed texture (|44| ) can also arise when the zeroes of the lopsided Dirac matrix 
in (RO) are replaced by small quantities. It is possible to construct models along this line 
based on a spontaneously broken U(1)f flavor symmetry, where 5, e and 7] are given by 
positive powers of one or more symmetry breaking parameters. Moreover, by playing with 
the U(l)p charges, we can adjust, to certain extent, the relative hierarchy between rj, e 
and 5 0, [43], [45], ||], |47|, |48| , as we will see in section 8. The texture ([44]) can also be 
generated in SUSY models with i?-parity violation |50| . 

After a first rotation by an angle 623 close to 7r/4 and a second rotation with #13 w e, 
we get 

e 0\ 



m 



( 5 + e 2 
e 




V 






2 



(45) 
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up to order one coefficients in the small entries. To obtain a large solar mixing angle, 
we need \rj — 8\ < e. In realistic models there is no reason for a cancellation between 
independent perturbations and thus we assume both 5 < e and 77 < e. 

Consider first the case 5 ~ e and 77 < e. The solar mixing angle 9 12 is large but not 
maximal, as preferred by the MSW LA solution. We also have Am^ (m ~ 4m 2 , Am 2 nra ~ 

Amltm^ 2 and 



AmLn ■ (46) 

If r] ps e and <5 e, we still have a large solar mixing angle and Am 2 un ps e 2 AraJ (m , as 
before. However m ee will be much smaller than the estimate in (|46|) . Unfortunately, this 
is the case of the models based on the above mentioned U(1)f flavor symmetry that, at 
least in its simplest realization, tends to predict 5 ~ e 2 . In this class of models we find 

/ Am 2 \ * 



V Am 



atm , 



below the sensitivity of the next generation of planned experiments. It is worth to mention 
that in both cases discussed above, we have 



7 13 



1 



/AmjLV 

yAm 2 im y 



(48) 



which might be very close to the present experimental limit. 

If both 5 and 77 are much smaller than e, the 12 block of m u has an approximate 
pseudo-Dirac structure and the angle 9 12 becomes maximal. This situation is typical of 
some models where leptons have U(1)f charges of both signs whereas the order param- 



eters of U(1)f breaking have all charges of the same sign [47]. We have two eigenvalues 
approximately given by ±m e. As an example, we consider the case where 77 = and 
5 m e 2 . We find sin 2 2# 12 ~ 1 — e 2 /4, Am 2 Mn m m 2 e 3 and 



In order to recover the MSW LA solution we would need a relatively large value of e. This 
is in general not acceptable because, on the one hand the presence of a large perturbation 
raises doubts about the consistency of the whole approach and, on the other hand, in 
existing models where all fermion sectors are related to each other, e is never larger than 
the Cabibbo angle. Therefore, the case 5, 77 <C e can be more easily adapted to fit the LOW 
solution where the solar frequency is small. As a consequence, m ee is beyond the reach of 
the next generation of experiments, whereas #13 might be tested at future facilities. 



5.4 Summary 

Given the present experimental knowledge, which favours Am 2 4m as the leading oscillation 
frequency and two large mixing angles, 6*23 and 6*i 2 , it is natural to define a zeroth order 
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approximation of the theory, where /S.m? sun and #13 vanish (which allows us to neglect 
the CP-breaking parameter (p) whereas #23 and 9 12 are maximal. For each pattern of 
neutrino masses, we have considered the most interesting textures that arise in this limit. 
This approximation is of course not realistic and should be regarded only as a limiting 

case, possibly arising from an underlying symmetry. Many effects can perturb this limit, 
such as small symmetry breaking terms, radiative corrections, effects coming from residual 
rotations needed to diagonalize the charged lepton mass matrix, or to render canonical the 
leptonic kinetic terms. Some of these will be discussed more in detail in the next sections. 
It turns out that in most of the existing models the leading textures are modified by small 
perturbations having a simple structure, such as those called Dl, D2, I and N in eqs. (|22|), 
(0), © and ©. 

We have analyzed these perturbations, with the hope that the results are sufficiently 
representative of the many existing models. Of course there is no guarantee that this 
discussion can cover all theoretical possibilities. Moreover, if /S.m 2 sun and #13 were as large 
as experimentally allowed, the perturbations would become large, the whole approach 
could become questionable and the data would be more appropriately described by an 
anarchical framework. 

A remarkable feature is that most models continue to predict an almost maximal so- 
lar angle, even after inclusion of the perturbations. This is often due to an approximate 
pseudo-Dirac structure in the 12 sector, which, at leading order, forces #12 = vr/4. Excep- 
tions to this trend are given by some of the possibilities offered by the normal hierarchy, 
for which 9\2 is undetermined at leading order. 

It is also apparent that, apart from the case of degenerate spectrum realized with a 
texture similar to the one of flavor democracy (D2), the possibility of measuring m ee with 
the next generation of experiments seems to be significant only if the solar oscillation 
frequency is very close to the upper part of the range allowed by the MSW LA solution. 



6 Importance of Neutrinoless Double Beta Decay 

The discovery of Qvj3j3 decay would be very important because it would establish lepton 
number violation and the Majorana nature of z/'s. Indeed oscillation experiments cannot 
distinguish between pure Dirac and Majorana neutrinos. Moreover, the search for Qvf3f3 
decay provides information about the absolute spectrum, while neutrino oscillations are 
only sensitive to mass differences. Complementary information on the sum of neutrino 
masses is also provided by the galaxy power spectrum combined with measurements of 



the cosmic microwave background anisotropics ||51|| . As already mentioned the present 
limit from Qv(3j3 is \m ee \ < 0.2 eV or to be more conservative \m ee \ < 0.3 — 0.5 eV [j3(i[ . 
In this respect it is interesting to see what is the level at which a signal can be expected 
or at least not excluded in the different classes of models [0, [53|] . For 3- neutrino models 
with degenerate, inverse hierarchy or normal hierarchy mass patterns, starting from the 
general formula in eq. (0), it is simple to derive the following bounds. 
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Table 2: Order-of-magnitude estimates for #13. 
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a) Degenerate case. If |m| is the common mass, apart from a phase, and taking S13 = 0, 
which, as already observed, is a safe approximation in this case, we have m ee = 
|m|(c 2 2 ± s\ 2 ). Here the phase ambiguity has been reduced to a sign ambiguity 
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which is sufficient for deriving bounds. So, depending on the sign we have m ee = \m\ 
or m ee = \m\cos26i2. We conclude that in this case m ee could be very close to the 
present experimental limit because \m\ can be sizeably larger than the Ov/3/3 bound 
(|m| < O(leV)), but should be at least of order 0(y / Am^ m ) ~ O(10~ 2 eV) unless 
the solar angle is practically maximal, in which case the minus sign option can be 
as small as required. 

b) Inverse hierarchy case. In this case the same approximate formula m ee = \m\(c^ 2 ± 
s\ 2 ) holds because is small and S13 can be neglected. The difference is that here 
we know that \m\ « \J Aml tm so that \m ee \ < \J Aml tm ~ 0.05 eV. 

c) Normal hierarchy case. Here we cannot in general neglect the m.3 term. However in 
this case \m ee \ ~ J Am^ m s\ 2 + \J Am^ m s\ z and we have the bound \m ee \ < a 
few 10~ 3 eV. 

Recently evidence for Ovflfl was claimed in ref. at the 2-3<r level with \m ee \ ~ 0.39 eV. 
If confirmed this would rule out cases b) and c) and point to case a) or to models with 
more than 3 neutrinos. 



7 Expectations for #13 

The measurement of #13 represents one of the main challenges for the next generations 
of experiments on neutrino oscillations, which, with the help of very intense neutrino 



beams [po] , might reach a sensitivity of few percent on #13. A sizeable #13 would have 
an important impact on the observability of CP-violating effects in the leptonic sector. 
We collect in table |2] our estimates of #13 for the various textures considered in section 5. 
Figure |I] displays the expectations for 6*13. Such expectations are of course very rough and 
are not meant to be statistically meaningful. It is however interesting to note that the 
MSW LA solution favours #13 in an experimentally accessible range, for all textures but for 
the inverted hierarchy, where we find a strong suppression. The LOW solution prefers a 
smaller, unobservable 6*13, with the possible exception of the texture D2, corresponding to 
flavour democracy. We recall that if the neutrino mass matrix is structureless as advocated 
by the anarchical framework, then # 13 is naturally expected to be very close to its present 
experimental bound. 



8 Grand Unified Models of Fermion Masses 

We have seen that the smallness of neutrino masses interpreted via the see-saw mechanism 
directly leads to a scale A for L non-conservation which is remarkably close to Mqjjt- 
Thus neutrino masses and mixings should find a natural context in a GUT treatment 
of all fermion masses. The hierarchical pattern of quark and lepton masses, within a 
generation and across generations, requires some dynamical suppression mechanism that 
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Figure 1: Order-of-magnitude estimates of #13, for MSW LA solution (red) and LOW 
solution (blu). There is no difference between MSW LA and LOW for the texture D2. 
The upper limits have been obtained by using the approximate expression of table 1 and 
the values: Am 2 atm = 2.5 ■ 1CT 3 eV 2 , Am 2 sun = 6.0 ■ 10" 5 eV 2 (10" 7 eV 2 ) for MSW LA 
(LOW). The continuous and dashed lines show the present upper bound on #13 and the 
possible reach of long baseline experiments with high intensity neutrino beams. 

acts differently on the various particles. This hierarchy can be generated by a number 
of operators of different dimensions suppressed by inverse powers of the cut-off A c of the 
theory. In some realizations, the different powers of 1/ A c correspond to different orders in 
some symmetry breaking parameter Vf arising from the spontaneous breaking of a flavour 
symmetry. In the next subsections we describe some simplest models based on SU(5) x 
U(l)p and on SO(10) which illustrate these possibilities. 

8.1 Models Based on Horizontal Abelian Charges 

We discuss here some explicit examples of grand unified models in the framework of a uni- 
fied SUSY SU(5) theory with an additional U(1)f flavour symmetry. The SU(5) generators 
act "vertically" inside one generation, while the U(1)f charges are different "horizontally" 
from one generation to the other. If, for a given interaction vertex, the U(l)p charges 
do not add to zero, the vertex is forbidden in the symmetric limit. But the symmetry is 
spontaneously broken by the VEV's Vf of a number of "flavon" fields with non vanishing 
charge. Then a forbidden coupling is rescued but is suppressed by powers of the small 
parameters f//A c with the exponents larger for larger charge mismatch f56f . We expect 
Vf ^ Mqut and A c ^ Mpi. Here we discuss some aspects of the description of fermion 
masses in this framework. 

In these models the known generations of quarks and leptons are contained in triplets 
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and , (i = 1, 2, 3) corresponding to the 3 generations, transforming as 10 and 5 of 
SU(5), respectively. Three more SU(5) singlets describe the RH neutrinos. In SUSY 
models we have two Higgs multiplets, which transform as 5 and 5 in the minimal model. 
The two Higgs multiplets may have the same or different charges. We can arrange the unit 
of charge in such a way that the Cabibbo angle, which we consider as the typical hierarchy 
parameter of fermion masses and mixings, is obtained when the suppression exponent is 
unity. Remember that the Cabibbo angle is not too small, A ~ 0.22 and that in U(1)f 
models all mass matrix elements are of the form of a power of a suppression factor times a 
number of order unity, so that only their order of suppression is defined. As a consequence, 
in practice, we can limit ourselves to integral charges in our units (for example, vA ~ 1/2 
is already almost unsuppressed) . 

There are many variants of these models: fermion charges can all be non negative with 
only negatively charged flavons, or there can be fermion charges of different signs with 
either flavons of both charges or only flavons of one charge. We can have that only the top 
quark mass is allowed in the symmetric limit, or that also other third generation fermion 
masses are allowed. The Higgs charges can be equal, in particular both vanishing or can 
be different. We can arrange that all the structure is in charged fermion masses while 
neutrinos are anarchical. 



8.1.1 F(fermions)> 



Consider, for example, a simple model with all charges of matter fields being non negative 
and containing one single flavon 9 of charge F= —1. For a maximum of simplicity we also 
assume that all the third generation masses are directly allowed in the symmetric limit. 
This is realized by taking vanishing charges for the Higgses and for the third generation 
components ^3°, and ^>\. For example, if we define F(\l/f) = qf 1 (R = 10,5, 1; i — 
1, 2, 3) we could take @ |§ |§, g7| (see also [0]) 



(g 10 ~ 10 



1 >?r>93°) 



(3,2,0) 
(2,0,0) 



(50) 



A generic mass matrix has the form 



m 



(yii^ qi+q[ yi2\ qi+q ' 2 yi 3 X 11+q * 

y 2 i\ q2+q 'i y22\ q2+q * y2^ q2+q 'z 1 r 
\y 3l X^ y 32 \ q3+q * y33* q3+q *. 



(51) 



10 \ 



where all the are of order 1 and (qi, q'j) are the charges of (\l/ 10 , \I/ 10 ) for m u , of (\l/ 5 , ^ 
for rrid or mf, of , \1/ 5 ) for (the Dirac v mass), and of (\E fl , ^Z 1 ) for M, the Majorana 
v c mass, respectively. We have A = ($)/ A c and the quantity v represents the appropriate 
VEV or mass parameter. It is important to observe that m can be written as: 



rn 



\ny\Riv 



(52) 
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where A# = Diag(X q ^ , X q ? , X q s) and y is the matrix. The models with all non nega- 
tive charges and one single flavon have particularly simple factorization properties. For 
example, if we start from the Dirac v matrix: mp = XiyoX^Vu and the v c Majorana 
matrix M = Xiy^XiA and write down the see-saw expression for m v = m^M -1 m£>, 
we find that the dependence on the q 1 charges drops out and only that from q 5 remains. 
As a consequence the effective light neutrino Majorana mass matrix m u can be written in 
terms of q 5 only: m v = X^(yJ ) yJ / lyD)X^ > -v^/A. In addition, for the neutrino mixing matrix 
Uij, which is determined by m u in the basis where the charged leptons are diagonal, one 
can prove that ~ A' 9i_9 J 5 ', in terms of the differences of the 5 charges, when terms that 
are down by powers of the small parameter A are neglected. Similarly the CKM matrix 
elements are approximately determined by only the 10 charges ]5(J: V£j KM ~ A' 9i '. 



With these results in mind, we understand that the q 10 charge assignments in eq. (|50|) are 
determined by requiring V us ~ A, V c b ~ A 2 and V u b ~ A 3 . However the same g 10 charges 
also fix m u : m c : rrit ~ A 6 : A 4 : 1. The experimental value of m u (the relevant mass values 
are those at the GUT scale: m = m{MGUT) @) would rather prefer q\° = 4. Taking 
into account this indication and the presence of the unknown coefficients y^ ~ 0(1) it is 
difficult to decide between q{° = 3 or 4 and both are acceptable. 

Turning to the 5 charges, the entries q\ = gf = have been selected in eq. (|50|) so 
that the 22, 23, 32, 33 entries of the effective light neutrino mass matrix m v are all 0(1) 
in order to accommodate the nearly maximal value of S23. The small non diagonal terms 
of the charged lepton mass matrix cannot change this. In fact, for qj°, qf chosen as in eqs. 
TUT) we obtain: 




m d = A 3 A 2 1 \v d = (rm) 1 , m u = A 2 1 1 ~T » (53) 




where v u ^ are the VEVs of the Higgs doublets. Note that the patterns : m s : nib 
m e : m M : m T ~ A 5 : A 2 : 1 are acceptable (but also q\ = 3 would be possible). One 
difficulty is that for m v the subdeterminant 23 is not suppressed in this case, so that the 
splitting between the 2 and 3 light neutrino masses is in general small. In spite of the fact 
that rriD is, in first approximation, of the form in eq. (|39|) the strong correlations between 
mo and M implied by the simple charge structure of the model destroy the vanishing of 
the 23 subdeterminant that would be guaranteed for generic M. Models of this sort have 
been proposed in the literature J58L f45|. The hierarchy between m 2 and m 3 is considered 



accidental and better be moderate. The preferred solar solution in this case is MSW SA 
because if mi is suppressed (and some suppression is needed if we want S13 small) the 
solar mixing angle is typically small. However, if with a moderate fine tuning we stretch 
by hand 777,2 to become sufficiently close to mi then the MSW LA solution could also be 
reproduced. From eq. (|53D, taking v u ~ 250 GeV, the mass scale A of the heavy Majorana 
neutrinos turns out to be close to the unification scale, A ~ 10 15 GeV. 

A different interesting possibility |6(| is to recover an anarchical picture of neutrinos 
by taking {q™, g 2 10 , g 3 10 ) = (4, 2, 0) and (qf' 1 , gf' 1 , qf 1 ) = (0, 0, 0). The 10 charges lead to an 
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acceptable pattern for the m u matrix and Vckmi as already discussed. For down quarks 
and charged leptons we obtain a weakened hierarchy, essentially the square root than that 
of up quarks: : m s : m& ~ m e : m^ : m T ~ A 4 : A 2 : 1. Finally in the neutrino sector 
the anarchical model is realized (both and M are structureless). 

Note that in all previous cases we could add a constant to qf, for example by taking 
(?i, Q21 ?I) = (4,2,2). This would only have the consequence to leave the top quark as 
the only unsuppressed mass and to decrease the resulting value of tan/3 = v u /vd down 
to \ 2 m t /m h . A constant shift of the charges q\ might also provide a suppression of the 
leading v c mass eigenvalue, from A c down to the appropriate scale A. One can also consider 
models where the 5 and 5 Higgs charges are different, as in the "realistic" SU(5) model of 
ref. |61|]. Also in these models the top mass could be the only one to be non vanishing in 
the symmetric limit and the value of tan (3 can be adjusted. 

8.1.2 F(fermions) and F(flavons) of both signs 

Models with naturally large 23 splittings are obtained if we allow negative charges and, 
at the same time, either introduce flavons of opposite charges or stipulate that matrix 
elements with overall negative charge are put to zero. For example, we can assign to the 
fermion fields the set of F charges given by: 

(ql°,ql°,ql°) = (3,2,0) 

(qiqii) = (6,0,0) 6>2a>o 

(ql,ql,ql) = (a,-a,o) . (54) 

We consider the Yukawa coupling allowed by U(l)p-neutral Higgs multiplets in the 5 and 
5 SU(5) representations and by a pair 9 and 6 of SU(5) singlets with F= 1 and F= — 1, 
respectively. If b = 2 or 3, the up, down and charged lepton sectors are not essentially 
different than in the previous case. Also in this case the 0(1) off-diagonal entry of mi, 
typical of lopsided models, gives rise to a large LH mixing in the 23 block which corresponds 
to a large RH mixing in the d mass matrix. In the neutrino sector, the Dirac and Majorana 
mass matrices are given by: 








(\ 2a 


1 


A a \ 




M = 




A ,2a 


A ,a A 






V A a 


A ,a 





m D = A b " a A ,a A ,a \v u , M = 1 \ lZa A ,a A , (55) 



where A' is given by (0)/A c and A as before denotes the large mass scale associated to 
the RH neutrinos: A ^> v u ^. After diagonalization of the charged lepton sector and after 
integrating out the heavy RH neutrinos we obtain the following neutrino mass matrix in 
the low-energy effective theory: 




m v = X b 1 + A a A ,a 1 + A a A ,a -f . (56) 
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The 0(1) elements in the 23 block are produced by combining the large LH mixing induced 
by the charged lepton sector and the large LH mixing in mp. A crucial property of m v is 
that, as a result of the see-saw mechanism and of the specific U(1)f charge assignment, the 
determinant of the 23 block is automatically of 0(X a X' a ) (for this the presence of negative 
charge values, leading to the presence of both A and A' is essential PS| , fT7|). The neutrino 
mass matrix of eq. ( [561) is a particular case of the more general pattern presented in eq. 
(f44[), for J « A 26 , e « A 6 and r] m X a X' a . If we take A rs A', it is easy to verify that the 
eigenvalues of m v satisfy the relations: 

mi : m 2 : m 3 = X 2{b ~ a) : X 2a : 1 . (57) 

The atmospheric neutrino oscillations require m 2 ~ 10~ 3 eV 2 . The squared mass difference 
between the lightest states is of 0(X 4a ) m|, not far from the MSW solution to the solar 
neutrino problem if we choose a — 1. In general Z/ e3 is non-vanishing, of 0(X b ). Finally, 
beyond the large mixing in the 23 sector, m u provides a mixing angle 9 12 ~ X b ~ 2a in the 12 
sector. For b > 2a, we recover a small solar mixing angle. For instance, taking 6 = 3 and 
a — 1, 9±2 becomes close to the range preferred by the MSW SA solution. When b = 2a, 
as for instance in the case b = 2 and a = 1, the MSW LA solution can be reproduced. 



8.1.3 F(fermions) of both signs and F(flavons) <0 

A general problem common to all models dealing with flavour is that of recovering the 
correct vacuum structure by minimizing the effective potential of the theory. It may be 
noticed that the presence of two multiplets 9 and 9 with opposite F charges could hardly 
be reconciled, without adding extra structure to the model, with a large common VEV for 
these fields, due to possible analytic terms of the kind {99) n in the superpotential. We find 
therefore instructive to explore the consequences of allowing only the negatively charged 
9 field in the theory. 

It can be immediately recognized that, while the quark mass matrices previously dis- 
cussed are unchanged, in the neutrino sector the Dirac and Majorana matrices are obtained 
from eq. (|55|) by setting A' = 





m D = X b ~ a \ v u , M = 1 A . (5 



The zeros are due to the analytic property of the superpotential that makes impossible 
to form the corresponding F invariant by using 9 alone. These zeros should not be taken 
literally, as they will be eventually filled by small terms coming, for instance, from the 
diagonalization of the charged lepton mass matrix and from the transformation that put 
the kinetic terms into canonical form. It is however interesting to work out, in first 
approximation, the case of exactly zero entries in and M, when forbidden by F. The 
neutrino mass matrix obtained via see-saw from mo and M has the same pattern as the 
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one displayed in eq. (|56|). A closer inspection reveals that the determinant of the 23 block 
is identically zero, independently from A. This leads to the following pattern of masses: 

mi : m 2 : m 3 = X b : X b : 1 , m \-m\ = 0(A 3fe ) . (59) 

Moreover the mixing in the 12 sector is almost maximal: 

#i2 = ^ + 0(A fe ) . (60) 

For 6 = 3 and A ~ 0.2, both the squared mass difference [m\ — m^)/m\ and sin 2 2# 12 
are remarkably close to the values required by the VO solution to the solar neutrino 
problem. This property remains reasonably stable against the perturbations induced by 
small terms (of order A 5 ) replacing the zeros, coming from the diagonalization of the 
charged lepton sector and by the transformations that render the kinetic terms canonical. 
By choosing b = 2 we obtain the LOW solution. We find quite interesting that also 
the just-so and the LOW solutions, requiring an intriguingly small mass difference and a 
bimaximal mixing, can be described, at least at the level of order of magnitudes, in the 
context of a "minimal" model of flavour compatible with supersymmetric SU(5). In this 
case the role played by supersymmetry is essential, a non-supersymmetric model with 9 
alone not being distinguishable from the version with both 9 and 9, as far as low-energy 
flavour properties are concerned. 

In conclusion, models based on SU(5) x U(1)f are clearly toy models that can only aim 
at a semiquantitative description of fermion masses. In fact only the order of magnitude 
of each matrix entry can be specified. However it is rather impressive that a reasonable 
description of fermion masses, now also including neutrino masses and mixings, can be 
obtained in this simple context, which is suggestive of a deeper relation between gauge 
and flavour quantum numbers. There are 12 mass eigenvalues and 6 mixing angles that 
are specified, modulo coefficients of order 1, in terms of a bunch of integer numbers (from 
half a dozen to a dozen), the charges, plus 1 or more scale parameters. In the neutrino 
sector we have seen that the scheme is flexible enough to accommodate all the solutions 
that are still possible. Models aiming at a realistic unification of electroweak and strong 
interactions should of course address other important questions such as the doublet-triplet 
splitting and its stability against quantum corrections, a proton lifetime compatible with 
the existing limits, a correct gauge coupling constant unification and the consistency with 
present bounds on flavour violation p2| . Encompassing all these features in a consistent 
and possibly simple model is a formidable task, that might require to go beyond the 
conventional formulation in terms of a four- dimensional quantum field theory. 

8.2 GUT Models based on SO(10) 

Models based on SO(10) times a flavour symmetry are more difficult to construct because a 
whole generation is contained in the 16, so that, for example for U(l)p, one would have the 
same value of the charge for all quarks and leptons of each generation, which is too rigid. 
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But the mechanism discussed so far, based on asymmetric mass matrices, can be embedded 

The 33 entries 



in an SO(10) grand-unified theory in a rather economic way |26|, 44, 33 1 
of the fermion mass matrices can be obtained through the coupling I63I63IO// among 
the fermions in the third generation, I63, and a Higgs tenplet 10h- The two independent 
VEVs of the tenplet v u and Vd give mass, respectively, to tjv T and b/r. The key point to 
obtain an asymmetric texture is the introduction of an operator of the kind 16216^16316^ 
. This operator is thought to arise by integrating out an heavy 10 that couples both to 
16 2 16# and to 16 3 16' H . If the 16# develops a VEV breaking SO(10) down to SU(5) at 
a large scale, then, in terms of SU(5) representations, we get an effective coupling of the 
kind 82 1038^ , with a coefficient that can be of order one. This coupling contributes to 
the 23 entry of the down quark mass matrix and to the 32 entry of the charged lepton 
mass matrix, realizing the desired asymmetry. To distinguish the lepton and quark sectors 
one can further introduce an operator of the form 16jl6jl0#45//, = 2,3), with the 
VEV of the 45h pointing in the B — L direction. Additional operators, still of the type 
16 j IQjlQn 16'^ can contribute to the matrix elements of the first generation. The mass 
matrices look like: 



m„ 



m D 
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a + e 
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Vd 



(61) 



(62) 



They provide a good fit of the available data in the quarks and the charged lepton sector 
in terms of 5 parameters (one of which is complex). In the neutrino sector one obtains a 
large 6*23 mixing angle, sin 2 29 12 ~ 6.6 • 1CT 3 eV 2 and 9 13 of the same order of 6*12. Mass 
squared differences are sensitive to the details of the Majorana mass matrix. 

Looking at models with three light neutrinos only, i.e. no sterile neutrinos, from a more 
general point of view, we stress that in the above models the atmospheric neutrino mixing 
is considered large, in the sense of being of order one in some zeroth order approximation. 
In other words it corresponds to off-diagonal matrix elements of the same order of the 
diagonal ones, although the mixing is not exactly maximal. The idea that all fermion 
mixings are small and induced by the observed smallness of the non diagonal Vckm matrix 
elements is then abandoned. An alternative is to argue that perhaps what appears to be 
large is not that large after all. The typical small parameter that appears in the mass 
matrices is A ~ \Jmd/m s ~ \fm^m T ~ 0.20 — 0.25. This small parameter is not so 
small that it cannot become large due to some peculiar accidental enhancement: either 
a coefficient of order 3, or an exponent of the mass ratio which is less than 1/2 (due for 
example to a suitable charge assignment), or the addition in phase of an angle from the 
diagonalization of charged leptons and an angle from neutrino mixing. One may like this 
strategy of producing a large mixing by stretching small ones if, for example, he/she likes 
symmetric mass matrices, as from left-right symmetry at the GUT scale. In left-right 
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symmetric models smallness of left mixings implies that also right-handed mixings are 
small, so that all mixings tend to be small. Clearly this set of models [35] tend to favour 
moderate hierarchies and a single maximal mixing, so that the MSW SA solution of solar 
neutrinos is preferred. 



9 Conclusion 

By now there are rather convincing experimental indications for neutrino oscillations. 
The direct implication of these findings is that neutrino masses are not all vanishing. 
As a consequence, the phenomenology of neutrino masses and mixings is brought to the 
forefront. This is a very interesting subject in many respects. It is a window on the 
physics of GUTs in that the extreme smallness of neutrino masses can only be explained 
in a natural way if lepton number conservation is violated. If so, neutrino masses are 
inversely proportional to the large scale where lepton number is violated. Also, the pattern 
of neutrino masses and mixings interpreted in a GUT framework can provide new clues on 
the long standing problem of understanding the origin of the hierarchical structure of quark 
and lepton mass matrices. Neutrino oscillations only determine differences of mf values 
and the actual scale of neutrino masses remain to be experimentally fixed. In particular, 
the scale of neutrino masses is important for cosmology as neutrinos are candidates for 
hot dark matter: nearly degenerate neutrinos with a common mass around 1- 2 eV would 
significantly contribute to Q m , the matter density in the universe in units of the critical 
density. The detection of Oi//3/3 decay would be extremely important for the determination 
of the overall scale of neutrino masses, the confirmation of their Majorana nature and 
the experimental clarification of the ordering of levels in the associated spectrum. The 
recent indication of a signal for 0i>(3(3 with m ee in a range around 0.4 eV, if confirmed, 
would point to a small but possibly non negligible contribution of neutrinos to Q m and, 
among models with 3 neutrinos, would favour those with a degenerate spectrum. The 
decay of heavy right-handed neutrinos with lepton number non-conservation can provide a 
viable and attractive model of baryogenesis through leptogenesis. The measured oscillation 
frequencies and mixings are remarkably consistent with this attractive possibility. 

While the existence of oscillations appears to be on a ground of increasing solidity, 
many important experimental challenges remain. For atmospheric neutrino oscillations 
the completion of the K2K experiment, now stopped by the accident that has seriously 
damaged the Superkamiokande detector, is important for a terrestrial confirmation of the 
effect and for an independent measurement of the associated parameters. In the near 
future the experimental study of atmospheric neutrinos will continue with long baseline 
measurements by MINOS, OPERA, ICARUS. For solar neutrinos it is not yet clear which 
of the solutions, MSW SA, MSW LA, LOW and VO, is true, although a preference for 
the MSW LA solution appears to be indicated by the present data. This issue will be 
presumably clarified in the near future by the continuation of SNO and the forthcoming 
data from KAMLAND and Borexino. Finally a clarification by MINIBOONE of the 



29 



issue of the LSND alleged signal is necessary, in order to know if 3 light neutrinos are 
sufficient or additional sterile neutrinos must be introduced, in spite of the apparent lack 
of independent evidence in the data for such sterile neutrinos and of the fact that attempts 
of constructing plausible and natural theoretical models have not led so far to compelling 
results. Further in the future there are projects for neutrino factories and/or superbeams 
aimed at precision measurements of the oscillation parameters and possibly the detection 
of CP violation effects in the neutrino sector. 

Pending the solution of the existing experimental ambiguities a large variety of the- 
oretical models of neutrino masses and mixings are still conceivable. Among 3-neutrino 
models we have described a variety of possibilities based on degenerate, inverted hierarchy 
and normal hierarchy type of spectra. Most models prefer one or the other of the possible 
experimental alternatives which are still open. It is interesting that the MSW LA solar 
oscillation solution, which at present appears somewhat favoured by the data, is perhaps 
the most constraining for theoretical models. For example, it is difficult to reproduce this 
solution in the inverted hierarchy models. The MSW LA solution can be obtained in the 
degenerate case, also including the anarchical scenario, and in the normal hierarchy case, 
but also in these cases rather special conditions must be met. In many cases the MSW 
LA solution corresponds to values of #13 rather large, not far from the present bound. The 
values of m ee (which determines the rate of 0v(3(3) that is found in models leading to the 
MSW LA solution are typically of order ^ Am^ n in the normal hierarchy case and can be 

even larger in the degenerate case, where the expected values are at least of order \J Am 2 atm 
(assuming that indeed for the MSW LA solution si 2 is large but not close to maximal). 

The fact that some neutrino mixing angles are large and even nearly maximal, while 
surprising at the start, was eventually found to be well compatible with a unified picture 
of quark and lepton masses within GUTs. The symmetry group at Mqut could be either 
(SUSY) SU(5) or SO(10) or a larger group. For example, we have presented a class of 
natural models where large right-handed mixings for quarks are transformed into large left- 
handed mixings for leptons by the transposition relation = m^T which is approximately 
realized in SU(5) models. In particular, we argued in favour of models with 3 widely 
split neutrinos. Reconciling large splittings with large mixing(s) requires some natural 
mechanism to implement a vanishing determinant condition. This can be obtained in 
the see-saw mechanism, for example, if one light right-handed neutrino is dominant, or 
a suitable texture of the Dirac matrix is imposed by an underlying symmetry. We have 
shown that these mechanisms can be naturally implemented by simple assignments of 
U(l)p horizontal charges that lead to a successful semiquantitative unified description of 
all quark and lepton masses in SUSY SU(5)x U(1)f- Alternative realizations based on 
the SO(10) unification group have also been discussed. 

In conclusion, the discovery of neutrino oscillations with frequencies that point to very 
small v masses has opened a window on the physics beyond the Standard Model at very 
large energy scales. The study of the neutrino mass and mixing matrices, which is still 
at the beginning, can lead to particularly exciting insights on the theory at large energies 
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possibly as large as M GUT . 
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